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Abstract. 

Time-series data from multiple gravitational wave (GW) detectors can be 
linearly combined to form a null-stream, in which all GW information will be 
cancelled out. This null-stream can be used to distinguish between actual GW 
triggers and spurious noise transients in a search for GW bursts using a network of 
{~f} ■ detectors. The biggest source of error in the null-stream analysis comes from the 

^J • fact that the detector data are not perfectly calibrated. In this paper, we present 

an implementation of the null-stream veto in the simplest network of two co- 
C ) ' located detectors. The detectors are assumed to have calibration uncertainties and 

\^J , correlated noise components. We estimate the effect of calibration uncertainties 

C ~ ) ■ in the null-stream veto analysis and propose a new formulation to overcome this. 

"*>~ ■ This new formulation is demonstrated by doing software injections in Gaussian 

^ , noise. 
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1. Introduction 



Given the time series data from a network of gravitational-wave (GW) detectors, one 
can find a particular linear combination of the data streams such that it does not 
contain any trace of GWs. The idea of this null-stream was proposed by Giirsel and 
Tinto in their classic work P^. Giirsel and Tinto proposed that the null-stream can 
be used to solve the 'inverse problem' of GW bursts, i.e., to compute the unknown 
quantities (two sky-positions and two polarizations) associated with the gravitational- 
waveform from the responses of three broad-band detectors. 

Recently, there has been a lot of interest in the null-stream among the GW 
community. The main reason for this rejuvenated interest is that the first generation 
of ground-based interferometric GW detectors El El El have started acquiring 
scientifically interesting data. Among the most promising astrophysical sources of 
GWs for these ground-based detectors are the transient, unmodelled astrophysical 
phenomena like supernovae explosions, Gamma-ray bursts and black hole/neutron 
star mergers - popularly known as 'unmodelled bursts'. Most of the algorithms 
currently used in burst-searches are time-frequency detection algorithms that look 
for short-lived excitations of power in the 'time-frequency map' constructed from the 
data [HI El IE] Since present-day interferometric GW detectors are highly complex 
instruments, the data often contains lots of noise transients which trigger the burst 
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detection algorithms. It is almost impossible to distinguish these spurious instrumental 
bursts from actual GW bursts using any physical model of the GW bursts. Thus, 
burst-data-analysis is usually performed as a coincidence analysis between multiple 
detectors. Although this 'coincidence requirement' considerably reduces the list of 
candidate burst triggers, one month of data can potentially produce hundreds of multi- 
detector random coincidences. So, it becomes absolutely necessary to have additional 
'waveform consistency tests' that distinguish actual GW bursts from spurious noise 
transients. A cross-correlation statistic that is formulated in |5] is already being used 
as a coherent waveform consistency test in the search for GW bursts in the data of 
LIGO detectors. This is making use of the fact that all the LIGO detectors are aligned 
(approximately) parallel to each other. 

Recently, it was proposed by [1()| that the null-stream can be used to distinguish 
between actual GW triggers and spurious noise transients in a search for GW bursts 
using any general network of detectors. The main idea is that, if the coincident 
triggers correspond to an actual GW burst, the null-stream constructed at the time 
of the triggers will contain no trace of the burst, and, will fall in to an expected noise- 
distribution. On the other hand, if the coincident triggers correspond to spurious 
instrumental bursts, the bursts will not necessarily cancel out in the null-stream, and 
the null-stream will contain some excess power. Many authors have proposed similar, 
but non-equivalent ways of implementing this. Our method is based on the excess 
power statistic UJ, which was first used by ^U] m this context. This method is formally 
developed and is studied in detail in |11| . For an alternative implementation, see |12) . 
A similar veto strategy using the null-stream constructed from the two calibrated 
output quadratures of GEO 600 0] detector is already being used to veto the burst 
triggers from GEO 600. This is discussed in [T3| , 

It was soon realised that the biggest source of error in the null-stream analysis 
comes from the fact that the detector data are not perfectly calibrated, due to various 
reasons. In such cases, the null-stream constructed from the data containing actual 
GW bursts will contain some residual signal and will deviate from the expected noise 
distribution. This paper tries to address such practical issues connected with the 
implementation of the null-stream veto in the burst-data-analysis using a network 
consisting of two co-located interferometric detectors, like the two LIGO detectors 
in Hanford, WA. The detectors are allowed to have calibration uncertainties and 
correlated noise components. In Section we briefly review the veto method in 
the case of two co-located detectors. In Sectional we estimate the effect of calibration 
uncertainties in the null-stream veto, and in Section |3 we lay out and demonstrate a 
formulation to overcome this effect. 

2. The null-stream veto 

In the case of detectors placed widely apart, the null-stream is a function of the 
antenna patterns, and hence, the source-position pp. But in the case of two co- located 
detectors, the null stream is particularly simple. If h\(t) and fi2(t) denote the properly 
calibrated time series data from the two detectors, the null-stream is just |10j : 

n{t) = hi(t) - h 2 (t), (1) 

or, in discrete notation 

rij = hij - h 2j , (2) 
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where hu and h^ are the sampled versions of hi(t) and hi{t). At the time of a set 
of coincident triggers, we construct the null-stream rij. Let Nf. denote the discrete 
Fourier transform (DFT) of rij computed using L samples of the data. We assume 
that the real and imaginary parts of Nk are derived from a multivariate Gaussian 
distribution of mean zero J and variance a\. Following jHIIHj, we compute the 'excess 
power' statistic from the null-stream: 

e= £ P ,, Pfe = ^iL. (3 ) 

7 ®~k 

k=m ^ 

It can be shown that e will follow a \ 2 distribution of 2M degrees of freedom in 
the case of a non-windowed DFT. But in the case of a windowed DFT, P^ are not 
independent \ 2 variables, and hence e will not follow a \ 2 distribution |14j . But,we 
note that the \ 2 distribution is a special case of the Gamma distribution. It can be 
shown that, to a very good approximation, e will follow a Gamma distribution with 
scale parameter a and shape parameter (3. These parameters are related to the mean 
[i t and variance of of the distribution of e by 

a=(^)\ P= a X (4) 

In order to estimate the parameters of the expected Gamma distribution, we generate 
a population of e from stationary data (i.e., data not containing the burst event under 
investigation, but surrounding it). To be explicit, we divide the data in to a number 
of segments each length L and compute e from each of these segments. From that 
population, ,u c and of can be estimated, and hence a and (3. 

It is known that the maximum signal-to-noise ratio (SNR) for the excess power 
statistic is achieved when the time-frequency volume used to compute the statistic 
is equal to the actual time- frequency volume of the signal [H]. Since the duration 
and bandwidth of the burst is estimated by the burst detection algorithm itself, this 
information is used to decide on the length (L) of the data used to compute Nk and 
the bandwidth over which Pk is summed over. 

If the e computed from the segment containing the burst is greater than a 
threshold, we veto the trigger. The threshold, t, giving a false-dismissal ('false- veto') 
probability of 7 can be found by 

/oo 
f(x;a,/3)dx, (5) 

where f(x; a, (3) is the probability density of the Gamma distribution with parameters 
a and (3. 

2.1. Software injections 

Let us define some terminology. The false-dismissal probability is the probability of 
an actual GW burst being falsely vetoed, and the 'false-veto fraction' is the fraction 
of GW bursts that are actually vetoed using this method. As a sanity check, we 
estimate the false-veto fraction by injecting some prototype gravitational-waveforms 
into two data streams of Gaussian white-noise and performing the analysis. If all of 

I If the real and imaginary parts of JVj. are derived from a non-zero-mcan Gaussian distribution, one 
can always convert them to mean-zero Gaussian variables by subtracting the sample mean /i^ from 
them. 
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our assumptions are true, the fraction of vetoed events among the injections should 
be equal to the chosen false-dismissal probability. 

The injected waveforms are Gaussian-modulated sinusoidal waveforms, of the 
form: 

A /2f 2 \ 1/4 
h(t) = h TS J^J shi[2iTf (t-t )}exp[-(t-t ) 2 /T 2 ] , (6) 

where Jo is the central frequency of the waveform (randomly chosen from the set {153, 
235, 361, 554, 849, 1053, 1245, 1534, 1856} Hz) and t is the time corresponding to 
the peak amplitude. We setup the envelope width asr = 2//o, which gives durations 
of approximately 1-20 ms. The corresponding quality factor is Q = \/2TrfoT = 8.9 and 
bandwidth is A/ = fo/Q ~ O.l/o- The quantity h rss is the root-sum-squared (RSS) 
amplitude: 

-1 1/2 



hr(t)&t 



= h ISS . (7) 



The h YSS is randomly chosen from the logarithmically-spaced interval (5 x 10 -22 , 1 x 
10~ 19 ). The amplitude spectral density (ASD) of the noise in the two data streams is 
chosen to be 1 x 10~ 22 /vT£z and 2 x 10~ 22 /\/Hz. We define the combined SNR, p, by 

P 2 = Pi + p\ , (8) 

where p\ and p2 are the optimal SNRs § in detecting the bursts in the two data 
streams, and quote this quantity while discussing the results. 

The veto analysis is performed with different thresholds. The fraction of vetoed 
events is plotted against the false-dismissal probability corresponding to the chosen 
threshold in Figure Q] (left). It can be seen that the estimated false-veto fraction is in 
very good agreement with the predicted false-dismissal probability. 

3. Calibration uncertainties 

So far, we have been assuming that the two data streams are perfectly calibrated. But, 
due to various limitations in the calibration procedure, the calibration of present-day 
interferometers can be subject to uncertainties of a few percent in amplitude and a 
few degrees in phase. This means that the null-stream constructed from the data 
containing actual GW triggers can contain some residual signal, and the computed 
statistic can vary from the expected distribution. This will result in a different false- 
dismissal probability than the one predicted by the hypothesis test. 

It can be shown that we see approximately the same residual signal power in the 
null-stream for both a 10% relative amplitude error and a 10 degree relative phase 
error; therefore, in principle, we need to consider both of these effects. However, it is 
not easy to conceive of a simple model for the possible relative phase error between two 
co-located detectors (since such errors most probably arise due to inaccuracies in the 
calibration process and will therefore be frequency dependent) ||. It is, however, easy 
to think of a simple model for one possible source of relative amplitude calibration 

§ The optimal SNR p in detecting a signal h(t) buried in the noise is defined by p 2 = 
4 In" \H(f)\ df/S n (f) where H(f) is the Fourier transform of the signal and S n {f) is the one-sided 
PSD of the detector noise. 

|| One another possible source of relative phase error would be a time-offset between the two data 
streams. Although it is unlikely that such an error exists at any significant level, it would be possible 
to include this in the following analysis if necessary. 
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Figure 1. [Left]: Estimated false- veto fraction plotted against the predicted false- 
dismissal probability, assuming that the two detectors are perfectly calibrated. 
[Right]: Estimated false- veto fraction in the presence of ±10% calibration 
uncertainty, plotted against the combined SNR of the injections, for three different 
values of the predicted false-dismissal probability (7). The three dashed lines show 
the predicted values of the false-dismissal probability in the absence of calibration 
uncertainties. 



error. Suppose we know the absolute calibration of each detector to some accuracy. 
Then we have the possibility that the difference between the two data streams contains 
a systematic calibration error that would be independent of frequency. We explore a 
way to deal with this type of error in the rest of this paper. 

As a simple model, we assume that the frequency-dependence of the calibration 
error is negligible, and that the calibration error is a constant scaling factor over short 
time-scales (of the order of seconds). In the context of this analysis, we can assume 
that one detector is perfectly calibrated and the other is calibrated with a wrong 
scaling factor, i.e., 



hi(t) =m(t)+h(t), 
h 2 (t)=n 2 (t) + (l + 6)h{t), 



(9) 



where n±(t) and n,2{t) are the detector noises in which a gravitational waveform h(t) 
is present. S is the relative calibration error which is assumed to be a real quantity 
and is constant over short time-scales. 

In order to estimate the effect of calibration error in the false-dismissal probability, 
injections are done simulating a relative calibration error of S = ±0.1 between the 
two data streams. The fraction of the vetoed events among the injections is shown 
in the right plot of Figure ^ as a function of the combined SNR of the injections. 
Different curves in the plot represent three different values of the predicted false- 
dismissal probability (7). It can be seen that the estimated false- veto fraction raises 
to alarmingly high values for strong signals. 



3.1. Dealing with calibration uncertainties 

In this section wc formulate a strategy to reduce the effect of calibration errors in the 
false-dismissal probability. We construct the following linear combination of the data 
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Figure 2. The excess-power statistic e'(£) constructed from null-stream n'(£;£) 
plotted against the free parameter £, for three different values of 8. We assume 
that the maximum possible absolute value of 5 is O.f. e' (—5/(1 + 5)) for each 
curve is marked with a square. It can be seen that the minimum value of e'(£) is 
less than, or equal to e'(— 5/(1 + 5)). 



streams by introducing a free parameter. £, in the null-stream construction, i.e., 

»'(t;0 = M*)-(l + 0>*a(*)- ( 10 ) 

Substituting for h\{t) and h2(t) from Eq.® gives, 

»'(t;0=n(*;0-(* + C + ^«)^(*), (11) 

where 

n(t;0="i(*)- (1 + "*(*). (12) 

is the 'perfect' null-stream. We can 'tunc' the parameter £ such that the residual 
signal disappears in Ea. (|llfl . This is accomplished by minimizing the 'excess power' 
in n'(t; £) by varying £ over an interval. 

As described in Section [2 the null-stream n(t; £) is divided into a number of 
short segments and the DFT of each segment is computed. It may be noted that 
in all segments except the one containing the burst, n'(t;£) = n(t;£), because the 
signal is absent in these segments. We denote the DFT of n(t; £) and n'(t; £) by Nk(Q 
and iVfc(^), respectively. The mean, //&(£), and variance, Ofc(£), of Nk(0 are estimated 
from the neighboring segments of the one containing the burst. These are the moments 
of the expected distribution of N' k (£) in the absence of the burst. Nk(£) is converted 
to a mean-zero Gaussian variable by subtracting the sample mean /Ufc(£) from it. The 
test statistic e'(£) is computed from the segment of n'(t; £) containing the burst: 

\K(0\ 2 



e'(0= J2 P ^,P k 

We minimize e'(£) by varying £ over an interval (£ m in,£max 
parameter-space can be fixed as 



Smin 



& 



(13) 
) ^f. The boundary of the 

(14) 



^f It is important to note that we minimise the excess-power, e'(£), in the null-stream, and not the 
total power, n' 2 (t;£). In the later case, if there are correlated noise in hi(t) and h,2(t), £ can take 
values which will minimise the correlated noise components in n'(t; £) instead of the residual signal. 
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were ±<5 max is the maximum allowed value of the calibration uncertainty. When 
£ — > —<5/(l + <5), the residual signal in the null-stream cancells out, and e'(£) — ► e(£), 
where 

m+M 



e(0 = L P <^) ' P * = T^T- ( 15 ) 

k=m a k^> 

which falls in to an expected Gamma distribution in the case of GW bursts. 

Since e is quadratic in N' k , apart from |iVfc| 2 and \(5 + £ + £6) Hk\ 2 , it contains 
also the cross-terms. This means that the minimum value of e'(£) could be less than 
e(£). This is illustrated in FigureEl for three different values of S. Thus, the obtained 
false-dismissal probability could be less than what is predicted by the hypothesis test. 
Since the cross-terms depend upon the actual value of (5, this adds an error-bar to 
the false-dismissal probability. But it may be noted that the actual false-dismissal 
probability is always less than (or equal to) what is predicted by the hypothesis test. 

3.2. Software injections with simulated calibration errors 

We generate two data streams with a simulated relative calibration error 5 according 
to Eq.JjJJ and perform the veto analysis, assuming that the |<5 max | = 0.1. The 
minimisation of e(£) is carried out using an optimised minimisation algorithm. The 
analysis is performed for three different values (-0.1, 0, 0.1) of 6 and the false-veto 
fractions corresponding to different thresholds arc estimated in each case (in all cases, 
we assumed that |# max | = 0.1). The mean value (of the three simulations) of the false- 
veto fraction corresponding to each threshold is plotted against the corresponding 
false-dismissal probability in Figure |3] (left). The extremum values corresponding 
to each threshold are used to generate the error-bars. It can be seen that the 
estimated false- veto fraction is always less than or equal to the predicted false-dismissal 
probability. 

The real figure-of-mcrit of a veto method is its ability to reject spurious events 
with a given false-dismissal probability. But, given that the probability density of the 
noise transients are not known a priori, there is no rigorous way of estimating the 
'rejection power' of the veto. The best we can do is to estimate the ability of the veto 
to reject a given glitch population. As a plausible estimation, we inject a population 
of sine-Gaussian waveforms with random parameters into two data streams. We then 
perform the veto analysis after choosing different thresholds. The estimated rejection 
power is plotted against the false-dismissal probability in Figure [3] (right). Since the 
'excess power' in the null-stream is proportional to the individual SNRs of the bursts in 
the two data streams, the rejection power is also proportional to the SNR. Each curve 
in the figure corresponds to a particular range of SNR p (see Eq.(JHJ). This suggests 
that veto efficiencies of > 90 % can be achieved with a false-dismissal probability of 
~ 1% for spurious noise transients with p > 10. 

4. Summary 

Th null-stream constructed from the data of multiple GW detectors can be used 
to distinguish between actual GW triggers and spurious noise transients in the 
search for GW bursts using a network of detectors. The biggest source of error 
in the analysis comes from the fact that the present-day detectors arc subject to 
calibration uncertainties. In this paper we have proposed an implementation of 
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Figure 3. [Left]: Estimated false- veto fraction plotted against the predicted 
false-dismissal probability. The relative calibration between hi(t) and h,2(t) is 
assumed to be between —0.1 < 5 < 0.1. [Right]: Rejection power of the veto 
plotted against the predicted false-dismissal probability. Each curve in the plot 
corresponds to a particular range of SNR. 



the null-stream veto in the search for GW bursts in the data of two co-located 
interferometers. We have estimated the effect of calibration uncertainties in the 
veto analysis by performing software injections in Gaussian noise with simulated 
calibration errors. A strategy is proposed to minimize this effect, assuming a simple 
model for the amplitude calibration-error and neglecting the errors in the phase 
calibration. This is done by introducing an additional free parameter in the null- 
stream combination and minimizing the excess-power in the null-stream. Finally, 
we compared the estimated fraction of falsely-vetoed GW-like injections with the 
predicted false-dismissal probability and found that the estimated fraction has a good 
agreement with the prediction. We also estimate the rejection power of the veto as a 
function of the false-dismissal probability. 
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